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Abstract 

Parton distributions contain factorizable final state interaction effects originating from the 
fast-moving struck quark interacting with the target spectators in deeply inelastic scattering. We 
show that these interactions give rise to gauge invariance of the transverse momentum-dependent 
parton distributions. As compared to previous analyses, our study demonstrates the existence of 
extra scaling contributions from transverse components of the gauge potential at the light-cone 
infinity. They form a transverse gauge link which is indispensable for restoration of the gauge 
invariance of parton distributions in the light-cone gauge where the gauge potential does not 
vanish asymptotically. Our finding helps to explain a number of features observed in a model 
calculation of structure functions in the light-cone gauge. 
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1 Parton model and QCD 



Hadron structure functions, measurable in deeply inelastic scattering, are genuine physical ob- 
servables which provide direct access to the microscopic constituents of matter and their intricate 
interaction dynamics. In the naive parton model the structure function is expressed in terms 
of a probability density q{x) to find a parton of a specific fiavor with a certain fraction x of the 
parent hadron's momentum. The underlying probabilistic picture for the scattering process relies 
on the fact that the constituents in a hadron boosted to the infinite momentum frame behave 
as a collections of noninteracting quanta due to time dilation. This simple and intuitive descrip- 
tion of hard reactions has found its firm foundation in rigorous field theoretical approach based on 
asymptotically free Quantum Chromodynamics (QCD). The result is factorization theorems which 
separate incoherent contributions responsible for physics of large and small distances involved in 
hard reactions in a universal and controllable manner: The physical observables such as structure 
functions are calculated as a convolution of QCD parton distributions in the hadrons and parton 
scattering cross sections. The parton model result arises as a lowest order term in the expansion 
in the coupling constant and inverse power of the hard momentum transfer of QCD factorization 
formulas. 

The QCD quark distribution follows from the factorization theorem in deeply inelastic scat- 
tering 

q{x) = \j {Pm^, 0)7+[e^, 0; 0, 0]^(0, 0)|P) , (1) 

where 

[e-, I; 0, i] ^ Pexp (^-ig di^A+{i^,i) 

is the gauge link between the quark fields, which arises from final state interactions between the 
struck quark and the target spectators. This interaction does not ruin factorization and is in fact 
much needed to maintain gauge invariance. On the other hand, the presence of this gauge link 
seems to spoil the interpretation of q{x) as a pure quark distribution, as the bilocal operator in 
the above expression is not obviously a quark number operator. The probabilistic interpretation 
is expected to hold only in the light-cone gauge [0, ||, 

A+ = 0, (2) 

since only the physical degrees of freedom remain with this choice. In this special gauge, the light- 
cone gauge link in the parton distribution in Eq. (|l]) formally vanish. The final state interaction 
effects between the outgoing parton and the target spectators have seemingly been removed and 
parton distributions depend only on the ground state wave functions. In particular, q{x) can be 
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calculated as the sum of the proton light-cone wave functions squared. The above observation is 
the basis for standard folklore that QCD parton distributions represent the density of partons in 
the target hadron or nucleus. 

Recently, however, the authors in Ref. questioned the validity of the standard result that the 
structure functions are given by the densities of partons in a target hadron. They argued that final 
state interactions between the struck quark and the target spectators have an intrinsic physical 
content and cannot be simply disposed off by a clever gauge choice. In the light-cone gauge, they 
found that although the final state interactions between the struck quark and the target spectators 
cancel after the transverse-momentum integrations, the effect remains in the interaction between 
the target spectators. They concluded on this ground that the parton distribution in the light- 
cone gauge is not given by Eq. (|l]): As such, structure functions cannot be just parton densities in 
an initial state hadron, which presumably does not know anything about final state interactions. 

The light-cone gauge analysis is complicated by the fact that the condition = alone 
does not remove all gauge freedom: x_-independent gauge transformations are still possible since 
they do not change the gauge condition. This leads to an difficulty in perturbation theory as the 
gauge propagator develops singularities in the Fourier conjugated momentum variable. In order 
to remove them, one has to fix the gauge completely. One way to do this is to adopt a specific 
boundary condition on components of the gauge potential. However, this may yield non-causal 
interactions, i.e., the final state interaction effects can be shifted entirely or partially either before 
or after the hard scattering. Moreover, the hadron wave functions do depend on the choice of the 
boundary conditions. The effect of the migration of final state interactions was noticed early in 
the analysis of shadowing 0, see also 0. 

Let us briefly summarize the results of this paper in light of findings in Ref. To fully 
exhibit the final state interaction effects, the Feynman parton distributions are considered as the 
integration of the transverse-momentum dependent parton distribution, 

q{x) = J (fk g(x, k) . 

As we already emphasized above, the final state interaction effects in q{x, k) arise from the inter- 
action between the active quark after it was hit by the virtual photon and the target remnant. 
We calculate it in the deeply inelastic lepton production of a jet with small transverse momentum 
and find 

q{x, k) = \f ^^e-«-+^'»-^(P|V^(e, [oo, oo; i-.i\^cl+ [oo, oo; 0, 0]c ^(0, 0) |P) , 
with path ordered exponentials stretched in light-like as well as transverse directions, 

[cx), oo;^_,^]c = [oo,oo;oo,^][oo,^;^_,^] . 



2 



This definition, which did not assume the vanishing of the vector potential at infinity, is completely 
gauge-invariant. In the light-cone gauge (0), the gauge link along the light-cone disappears, 
[oo,^;^_,^] = 1, and we are left with the link in the transverse direction, [oo, c»; oo, which 
is crucial to maintain the gauge independence under residual gauge transformations. The effect 
of final state interactions can thus be present even in the light-cone gauge. As is clear from the 
above equation, the usual transverse momentum-dependent parton distributions with light-cone 
gauge links only are incomplete when they are used in the light-cone gauge, — they are no longer 
gauge invariant under different choices of the gauge boundary conditions. The transverse link 
in our new definition is responsible for the final state interactions which generates the Sivers 
distribution function: the parton transverse momentum distributions asymmetry in a transversely 
polarized nucleon . This distribution is responsible fully for the single transverse-spin asymmetry 
discovered in and explained in [0. The mechanism which led to nonvanishing asymmetry 
was thought to be time-reversal odd and therefore prohibited in QCD. Our analysis supports the 
conclusion in that time-reversal does not impose any constraint on the Sivers function and 
hence the rich phenomenology of transverse momentum-dependent parton distributions as related 



to spin physics, see, e.g., 12, 13, 141 nB, 16, 17, 18 



With different choices of boundary conditions one can shift residual final state interactions from 
the final to the initial state. For the advanced boundary condition under which the transverse field 
A vanishes at = oo, the transverse gauge link vanishes as well, and the transverse momentum 
distributions are computable solely as a sum of the light-cone wave functions squared! Where are 
final state interactions now? They are now contained in the light-cone wave function due to the 
boundary condition on the gauge potential. Since the hadron wave function is a gauge-dependent 
object, it may acquire certain effects with the change of the gauge adopted to solve the bound 



state equation |T^. Choosing the light-cone gauge with a specific boundary condition one mimics 
final state interactions by a phase if in the wave function, = lale**^, which reproduces the effects 
of the gauge links if a different boundary condition is to be chosen. Therefore, the light-cone wave 
functions are special in that they encode not just the ground state properties of a system but also 
the final state interaction phases introduced through the choice of gauge conditions. This holds 
even though they represent stable bound state particles. This is not entirely surprising because 
final state interactions represented by the gauge links are quite general in nature, independent 
of details of the jet. In our view, this is a feature, rather than a defect, of the light-cone wave 
functions. 

Now return back to the ordinary Feynman distributions. The transverse link at infinity cancels 
when integration over k is performed. Therefore, independent of the residual gauge fixing, Eq. 
(1) is the correct formula for the parton distribution in any gauge. In the light-cone gauge, the 
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light-cone links reduces to unity, and the parton distribution can be calculated as the sum over 
the amplitude squared of the light-cone wave functions. The effects of final state interactions are 
partially cancelled through the transverse momentum integration and partially included already 
in the light-cone wave functions. So the parton distributions will be uniquely determined by the 
light-cone wave functions, q{x) ~ | and may be interepreted as the parton density in the 
spirit of the conventional Feynman parton model. We emphasize again here that the conclusion 
is independent of the residual gauge fixing in the light-cone gauge. 

Therefore the parton distributions do contain factorizable final state interactions. However, 
they can be completely hidden in the hadron wave functions if one chooses light-cone gauge and 
the advanced boundary condition. Then the parton distributions become densities in the sense 
that they can be calculated as the sum over the light-cone wave functions squared. As an example 
that the light-cone wave functions contain final-state intereaction effects, the light-cone amplitudes 
are complex when the advanced boundary condition is chosen [2^], and they do not possess simple 
time-reversal symmetry properties. Hence some single spin asymmetries can be attributed to, 
albeit in a non-causal way, the properties of the bound-state wave functions! 

A complementary view of the gauge links is to regard the (outgoing/incoming) "physical par- 
ton" as a free parton field with an attached gauge link extending from the parton position to 
(positive/negative) infinity. This "physical parton" field is gauge invariant and contains all the 
initial/final state interaction effects. The parton distributions are then the densities of these 
"physical partons" in a bound state. 

The outline of our presentation is as follows. In the next section we introduce objects of study 
and spell out our conventions. Next we address the issue of the gauge invariance of transverse 
momentum-dependent parton distribution in deeply inelastic scattering (DIS) by deriving the 
gauge link. As we stated above, it extends along the light-cone direction and is accompanied 
by an extra scaling contribution from the transverse components of the gauge field at light-cone 
infinity. The latter contribution would not be present for local gauge potentials, however, it affects 
the distribution in the light-cone gauge which generates non-local and non-causal interactions. 
In section ^, we calculate parton distributions in a simple scalar QED model of Ref. and 
demonstrate the effect of the final state interactions in the interplay with boundary conditions 
imposed on the gauge potential. We demonstrate the utter equivalence of different prescription 
for the light-cone pole on the amplitude level provided the transverse gauge link is accounted 
for. We show that with the advanced boundary conditions the final state interactions vanish and 
the parton distribution are defined solely by hadronic wave functions. In section ^ we show the 
equivalence of the dipole picture of deeply inelastic scattering in the aligned-jet kinematics, as 
discussed in Ref. 0, to the twist-two structure functions. In section ^ we address transverse 
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momentum- dependent parton distribution functions in the Drell-Yan (DY) process and discuss 
violation of the naive universahty with DIS distribution functions. A procedure for the complete 
fixing of the residual gauge freedom in the light-cone gauge @ by imposing a set of boundary 
conditions is given in the appendix 0. 



2 Definitions 

The main object of our study will be the structure functions which parametrize the hadronic part 
of the deeply inelastic scattering cross section and expressed, in turn, by the absorptive part of 
the forward Compton scattering amplitude 



W,, = ^'^mT,,. (3) 



The latter is a matrix element of the chronological product of quark electromagnetic currents 

T,. = ij d'ie^^<{P\T{j,{i)Um \P) . (4) 

Using the completeness of hadronic states Z]Ar|^)(^| = with a state |A^) consisting of 
particles with the total momentum P/v = J2k=i Pk ^i-nd certain quantum numbers, one gets 

X (P|j^(0)bi, ■ ■ ■,Pn){pi, ■ ■ ■,Pn\UO)\P) , (5) 

where — m^) = 6{po)S{p'^ — m?) imposes the on- mass-shell condition for an outgoing particle 
and the summation over N involves the summation over the number of particles populated the 
final states as well as their quantum numbers. 

For an unpolarized target, the hadronic tensor admits the following decomposition into the 
structure functions 



which depend on the photon virtuality = —q^ and the Bjorken variable x-q = Q"^/ {2q ■ P). In 
the Bjorken limit Q"^ oo with xb = fixed, these functions are expressed, 

F1{XB,Q') = ^F2{XB,Q') = IJ:QI (/.(^B,g') + fq{xB,Q')) , (7) 

in terms of QCD quark fg and antiquark fg distribution functions, at leading order of QCD 
perturbation theory and neglecting power-suppressed contributions. The dependence of fq^q on 
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arises only via the renormalization of the underlying composite operators and is logarithmic. 
In our subsequent considerations this will be irrelevant and, in order to simplify formulas, we 
neglect it in parton distributions. We also drop the summation over the quark flavors and set the 
quark charge to unity, Qq = 1. 

The Bjorken kinematics sets the microscopic dynamics on the light-cone, since the photon 
four-vector has large g_ = Q'^/{2x-b) and small g+ = — xb component so that the integrand in Eq. 
(H) does not oscillate provided ^+ ~ is small and ^_ ~ is large. In a frame where q does 
not have the transverse components, the light-cone decomposition of the external momenta reads 



We introduced the light-like vectors and n* such that 

= = , n ■ n* = 1 . (9) 

For an arbitrary four- vectors f^, we define 

= w • n , w_ = f • n* , V ■ u = v^U- + V-U+ + v± ■ u± . (10) 

We use Euclidean notations for transverse two-dimensional space with the metric Sa/s = —Qap = 
-(5'a/3 - naul - n*jaii) = diag(l, 1), so that 

v± ■ U± = —V ■ u . (11) 

Using these definitions one can easily project the hadronic structure functions via 



Q 



i3 



^1 = -7;9t.W,u , Fl = ^^W.. , (12) 



where we defined the function = F2 — 2xbFi. These structure functions are related to the 
absorption cross sections of transversely and longitudinally polarized photon by the hadron, Fi ~ 
(Tj. and Fl ^ a1 , respectively. In the Bjorken limit, Fi ^ Fl for spin-1/2 constituents, while for 
scalar partons the inequality is reversed, Fi <^ Fl- 



3 Parton distributions 

In this section we demonstrate how a gauge invariant parton densities arise in QCD description 
of the structure functions at leading twist. In particularly, we consider transverse momentum- 
dependent distributions. In QCD, the parton distributions are defined as hadronic matrix elements 
of quark bilocal operators. Since the parton fields enter at distinct space-time points, a gauge link 
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is needed to make the operators gauge invariant. This gauge hnk is generated in a hard scattering 
process by the final state interactions between the struck parton and the target remnants. Since 
the struck parton moves with a high energy, its interaction can be approximated through an 
eikonal phase. We show that the conventional light-cone link is not the only scaling contribution. 

A classic analysis demonstrated that non-gauge invariant quark bilocal correlators rep- 
resenting parton distributions in the zero-order approximation acquire path-ordered gauge links 
stretched along a non-light-like direction, nNL, from the position of quarks to infinity, and restore 
in this way their gauge invariance. A naive lifting of the direction to a light-like tinl = n was 
shown to lead to severe divergences [Q, which cancel only for integrated densities. Although, 
not implied by Ref. [Q, the following definition of the transverse- momentum dependent quark 
distributions has been unequivocally accepted in the literature 

where the path-ordered gauge link extends along the light-cone 

[oo,|;e-,|] = Pexp(^-2^7^"rfe-^+(^,0) • (14) 

Since the coordinate ^+ = 0, we do not display it, for brevity, here and the following presentation. 
As we will see, Eq. (p!3D is true only in a class of gauges where the gluon potential vanishes at 
^_ = 00. Note that the integrated parton distribution is just 

g(x) = J(fkq{x,k), (15) 

and generates quark and antiquark distributions q{x) = fq{x)9{x) — fqO{—x). 

For the photon scattering on a single parton with momentum £, we have for the hadronic 
tensor 

= ^E/ ^(27r)5+(Pj)(27r)^5(^) {Pn+Pj-P- q) {P\3,mPJ,N){pj,N\m\P) , 

(16) 

where p,j is the momentum of the observed final state quark (jet) after multiple rescattering with 
spectators N in the target fragments. The tree scattering amplitude corresponding to Fig. |1| (a) 
reads 

(pj,iV|j,(0)|P)(o) = u{pjh,{N\i^mP) , (17) 

where u is the Dirac spinor of the scattered quark. Substituting Eq. ([T7|) into (|T6D we get, after 
projection with (p!2D, the structure function Fi(xb,<5^) in the parton model 

i^i(a;B,Q') = \j (fpjq{x^.Pj) , (18) 

7 
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Figure 1: Multi-gluon attachments to the struck quark in deeply inelastic scattering which factorize 
to form the path-ordered exponential. 



in terms of the parton distribution 

1 fd^_ 



2 7 27r (27r)2 

As it stands the above correlation is not gauge invariant 



(19) 



3.1 Light-cone gauge link 

Consider now the contribution from the diagram in Fig. |I| (b). The one-gluon amplitude reads 



{Pj,N\j,{0)\P)^,)=gu{pj) 



where the free quark propagator is 
Momentum conservation gives 



Soik) 



-{N\4{h)So{pj-hhMO)\P) 



(20) 



(21) 



Pj = q + i + h 



so that the on-mass-shell condition Pj = results into 

£+ + ki+ = xb , 

in the Bjorken limit g_ —>■ oo. By expanding all vectors in Sudakov components (PH]), one can 
easily keep track of leading contributions in the scaling limit. For instance, the struck quark 
propagator contains a scaling term, 

1 7+ 



So{Pj - h] 



(1-1+ 



2g_(£+ -Xb) +iO 2ki+-i0 
with the remainder naively suppressed by extra powers of 1/q-- 



(22) 



8 



Making the light-cone decomposition of the gluon field 

= n*^A+ + n^A. + A^ , (23) 

one notices that the leading twist contribution comes from the first term on the right-hand side 
only, since the second one vanishes due to nilpotence of 7^ = 0, while the transverse component 
of the gauge field is of twist-three. Since the leading component in the quark density matrix u®u 
comes from the large term g_7+ we are allowed to replace 'U7_7+ by 2u on the right-hand side of 



Eq. (20). Finally, integrating with respect to ki we get 

/ Si;:^-^**') = ' L "^-^KM-o^o) . (24) 

Using these results, we find that the amplitude in one-gluon approximation differs from the three- 
level result (|l^ only by an extra factor, namely, 

{pj,N\UO)\P)ii) = u{PjhAN\i-w) / rfe~^+(e-,0,O)^(0)|P) . (25) 

<^ 

These considerations are generalizable in a straightforward manner to an arbitrary order. 
Namely, for n-gluon exchange, see Fig. |l] (c), we have 

{pj,N\jumP)in)=9MPj) J l[j^,W{k,)So{pj-k,)4{k,)So{pj-k^-k2)--- 



x4ikn)So \pj-Y.kj] 7.^(0) |P) . (26) 



Making repeatedly the same set of approximations as above, one simplifies the Dirac structure to 

" d% 1 1 1 



/II, 
n 



\ (27r)4 ki+ - iO ki+ + k2+-i0"' E"=i kj+ - tO 

xu{pj)j,{N\A+{k,)A+{k2) ■ ■ ■ A+{knMO)\P) , (27) 

where we have used the momentum conservation + X^ILi — ^b- The calculation of the 
momentum integrals is trivial and gives, similarly to Eq. (|24D, 



" d'^k- 11 1 

n (2^ fc^T^-^ fc^TTfc^T^-^ ■ ■ ■ T^:;k-^'^^^^'^'^^^^'^ ■ ■ ■ "^^^^'^^ 

/oo " 
n d^^+e{^^- - e(.+l)-)A+(ei-)A+(6-) • ■ ■ A+(a-) , (28) 

where = 0. Thus, 

_ _ POO Z"?!- 

(pj,iV|j,(0)|P)(„) = u{pj)^,{N\{-zgr / d^i- / d^2^--- 

Jo Jo 

X f^'""'' rfen-A+(ei-)^+(6-) ■ ■ • A+(a-)^(o)|p) . (29) 

Jo 
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Here, one immediately recognizes the nth term in the expansion of the path-ordered exponential 
Therefore, for the amphtude resummed to all orders, one gets 

(pj,iV|j,(0)|P) = ^(pj,iV|j,(0)|P)(„)=w(pj)7.(iV|Pexp / d^^A+{^,)) m\P) ■ (30) 

n=0 V "'0 / 

Substitution of this result into the hadronic tensor yields, indeed, the conventional quark distri- 
bution (1131). 



3.2 Transverse gauge link 

In the previous section we were not quite accurate and actually omitted contributions which scale 
in the Bjorken region. These terms survive only at a point of the momentum space and are, 
normally, assumed to be vanishing. However, this is not the case for all gauge potentials. And a 
counter-example is provided by the light-cone gauge where the potential is singular at the very 
same point. 

In the present circumstances it is more convenient to work in a frame similar to the Drell-Yan 
frame as was done previously in Refs. 0. Namely, the four-momentum of the current jet pj 



is light-like, p'j = 0, and can be chosen as one of the light-like vectors with the other one still fixed 
by the hadron momentum = n*^- Thus, we define 

n, ^ . (31) 

Obviously, = n*"^ = and n ■ n = 1. We decompose all Lorentz tensors via light-cone and 
transverse components defined in this basis 

Vf, = h^v^ + nlv+ + vj; , (32) 

and keep the same -|- index for contractions with v+ = n ■ v. In the Bjorken limit g_ oo, 
the difference between pj_ and g_ is negligible and, therefore, both frames, used before and here, 
coincide. 

Let us discuss first one-gluon exchange contribution in Eq. (^) . By looking at the denominator 
of the quark propagator one immediately notices that the scaling contribution in the Bjorken limit 
Pj_ oo, 

, ' , . (33) 

arises not only when one extracts the large pj_ component from the numerator but also when 
A;i+ ~ 1/pj- and keeps finite contributions in the numerator. A simple algebra gives 

u{pj)l,{^j - h)lu{N\Mh)m\P) ^ 2pj_u{pj)^,{N\A^{kMQ)\P) 

- u{pj)-i^}^,^,{N\A^{h)m\P) ■ (34) 
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where we have used the fact that, when keeping the transverse part of the gluon field, can 
now be pushed through giving zero when acting on the on-shell spinor due to the equation of 
motion u{j)j)'^j = 0. Thus, we find 

u{pj){N\4{h)Soipj-h)^MO)\P) ^ ^ .M PjhAN\A^ih)m\P) 



2pj_ki+ + kl - iO 



+ ^iPj) .^^ , ,2 _.j N\A^ih)m\P) ■ (35) 



The first term on the right-hand side is a contribution to the conventional light-cone link. When 
summed with contribution of multi-gluon exchanges it results as before into Eq. (pO]). We will drop, 
therefore, these terms completely and concentrate on effects of the second kind from transverse 
components of the gluon field. 

As we already emphasized before, in the scaling limit of pj_ — > oo, a finite contribution comes 
from = 0, i.e., when the exchanged gluon carries no longitudinal momentum. To see this, we 
exponentiate the denominator via the Chisholm representation 

1 

i 



/ (;Ae-'^(2p,/-fci++fc?-io)_ (3g) 
Jo 



Substituting it into (^) and then into Eq. (j20|), we can perform the integrations with respect to 
ki+ and which merely yield the Fourier transform of the gauge potential in these momentum 
components v4^(^- = l/(2Apj_),,^+ = 0, fci). Thus, in the scaling limit pj_ oo, the argument 
of v4(^) is set to = oo. The integration over A can now be trivially performed giving the 
propagator in the transverse space, 

{pj,N\MO)\P)ii) = 9u{pj) I §^^c.^^lAN\A^i^- = oo,e+ = 0,k,)m\P) ■ (37) 
To proceed further, we note that = oo, ^+ = 0, ^) must be a pure gaugeQ 

= oo, e+ = 0, 1) = v^m , (ss) 

since the field strength vanishes. The Fourier transform of this potential to the mixed, light-cone 
coordinate-transverse momentum representation, is Aq(^_ = oo,^+ = 0, A;) = ika4>{k), with 
being the Fourier transform of (f){^)- Substituting it into Eq. (^) we cancel the denominator k^, 
making use of 7^7/3 + 7/37a = -25a/3, and get 

(pj,iV|jUO)|P)(i) = -igu{pjhAPJ,N\mm\P) ■ (39) 



^Generally, a pure gauge field is A^{x) = {x)d^g{x) with a group-valued function g{x). We assume that g{x) 
represents small gauge transformations which are contractable to the unit clement g « 1 + 0. In the Eq. (|3^) we 
kept only the leading term in the perturbative expansion. 
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It is easy to see that 0(0) can be represented as a line integral 

POO 

0(0) = -/ d^-A{^,^). (40) 



The above expression is the first term in the expansion of an additional eikonal phase to the 
conventional one (|25|) , which has been neglected in the literature. Therefore, even in the light- 
cone gauge y4+ = 0, Fig. (b) does generate a non-zero contribution to the parton distribution. 

The momentum space procedure just outlined cannot be easily extended beyond single-gluon 
exchange. Therefore, the strategy will be to transform all factors in the integrand of Eq. (^) into 
the transverse coordinate space and do all manipulations there. To this end, we use 

^ = / rf2^e"=-% = -7^ / rf'^e^'^-^ V^lnl^l , (41) 



where we do not display a mass parameter which makes the argument of the logarithm dimen- 
sionless. Thus, the right-hand side of Eq. ( P7D can be equivalently written as 

-^^«(pj)7.7/37./^V^ln|^|(iV|V„0(O^(O)|P). (42) 

Next, one integrates by parts so that both derivatives act on the logarithm, Vq,V/3 In |^|. This is a 
symmetric tensor in its two-dimensional Lorentz indices so that the contracted Dirac matrices can 
be also symmetrized and reduced to the Kronecker symbol b^p by means of the Cliford algebra 
relation. Finally, since ln|^| is a two-dimensional Green function, one has V^ln|^| = 27r5''^^(^). 
These manipulations lead to the result in Eq. (39). 

For ra-gluon exchanges, we have the following contribution from the transverse components of 
the gluon field, similarly to Eq. (pT]), 

(p.,iVb.(0)|P)(„) = g-u{vjm J ^4(oo,fc,)^A^/04(oo,fe,) + 



(27r)2^ ^ ' kl -zO J {27ry^ ' ' (fci + fca)' - ^0 

iO 



^4((x,.fc„i ^-'7 — i.m\p)- («) 



The Fourier transformation gives 

{pj,N\jM\P)in) = Hg)MPj){N\J U^M^'^i)yM^i - ^2\Aioo,^,)y,\n\^, - ^,\. . . 

x4(oo, ^„_i)yn-i In - ^ J4(oo, In I^J 7um\P) ■ (44) 

Now, exploiting Eq. (|38|), we integrate by parts starting with and using y^i^ln|^^ ~ C2I = 
— In — ^2! = —2tt5^'^\^i — ^2)) then with respect to fi^st noticing that 

mAaiO = \^a<P'{0 , (45) 
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Figure 2: Graphical representation of transverse momentum-dependent parton distribution in the 
coordinate space with gauge hnks endowing nonlocal product of quark fields with gauge invariance. 

and then performing the same partial integrations as for ^j^. In this way one can integrate out all 
^'s. At the end, one gets 



As a final step, one simply notices that 



Jo Jo 



so that it forms a gauge link once resummed to all orders 



[00, 00; 00, 0] = P exp [W J ' A{oo, ^] 



(46) 



(47) 



(4J 



Therefore, restoring the light-cone gauge link, one finds the complete result for the amplitude 



(49) 



where 

[oo,oo;^_,^]c = [00, 00; oo,^][cx),^;^_,^] , (50) 

is shown in Fig. 0. 

Multiplying (^9]) by its complex conjugate, we deduce the gauge invariant transverse mo- 
mentum-dependent parton distribution 

1 fd^. d^$ 



q{x,k) 



je-'^^-^^'^-HPm, I) [00, 00; e-, 7+[oo, 00; 0, 0]c ^-(0, 0) |P) . (51) 



2J2n (27r)2 

The unitarity implies a partial cancellation of links at light-cone infinity 

[00, 00; ^]l[oo, 00; 0, 0]c = [00, ^; |]^[cx), I; 00, 0][oo, 0; 0, 0] 



(52) 



so that the definition (|T3D of the parton distribution, accepted in the literature, acquires an 
additional transverse link. 
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4 Parton distribution in scalar QED 



In this section we demonstrate, with a model calculation of parton distribution functions in scalar 
QED the role played by the transverse gauge link in the restoration of a gauge invariant 

answer on the amphtude level. In Ref. 0, this model was used to demonstrate the presence of 
non-pinched final state interactions in scaling contribution to the deeply inelastic cross section. We 
reformulate this in the language of parton distributions and show how the final state interaction 
can be shifted from the final to the initial state depending on the prescription on the spurious 
pole in the gluon propagator in the light-cone gauge. 

The model considered in Ref. [^] consists of heavy D and light charged scalars with masses 
M and m, respectively, interacting with massive U(l) gauge fields A^, 

/:sQED = + iT^.D)^ T^.D - - M^D^D - ^F'^, + , (53) 

via the covariant derivative = 9^ -|- igA^. Using the electromagnetic current of scalar QED 

J^ = i(p^ {d^,-d,}j(p, (54) 

in the hadronic tensor (p!6D, a simple calculation gives for the longitudinal structure function 
in the handbag approximation^] 

F^^-^XB, Q') = 2xB {n{xB) - it{-Xb)} , (55) 
with the scalar 0-parton distribution in the D-meson 

1 f ^e^-«-(p|<^t(o,o)[oo,oo;0,0][,za+[oo,oo;e-,0]c</>(e-,0)|P). (56) 
Z J Ztc 



TTiX] 



Here the subscript C denotes the path shown in Fig. |^. In practical computations we insert a 
complete set of intermediate states, so that the parton distribution reads 

vr(x) = y |^e-«~ X: / n ^(2^)^+ {pI - <) A(0)^iv(e-) , (57) 

where we introduced the amplitude 

An{^-) = (pi, . . . ,p;v| [oo, oo; 0]c <j){^-, 0)\P) . (58) 

In the following, we consider two-particle final state N = 2 with momenta pi = —p and 
P2 = P'. The amplitude expanded in perturbation series takes the form 

-^^(^-^ = / 7|^ie^'-^'"'^n2vr)^5(^)(^ + P'-P)T. 9'''A%\k) , (59) 



n=l 



^Recall that Fi ^ ax is zero for scalar quarks to this order. 
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Figure 3: One-gluon exchange contribution in the model calculation of scalar parton distribution 
in the Feynman gauge. 

where k is the total t-channel momentum 



n 



k = Y.ki, (60) 

i=l 

transferred from the target. For the small-x behavior, we are presently discussing, gluon t-channel 
intermediate states dominate, so that ki is the momentum of the ith gluon. Integrating out P^, 
k± and p+, we get for the small-x distribution function (|56|) 

2 



X f dp^ f d k d p 2 
The integrated components of momenta are set to 



E/("-i)aS;)(/c) 



n=l 



(61) 



B(p)_ ,^^E±^, (62) 

2p_ 2p_ 

by means of momentum conservation and on-mass-shell conditions. Here we introduced the con- 
vention 

D{p) = 2xp_ + p2 + - iO . (63) 

Now we are in a position to perform the computation of the amplitude in different gauges 
in order to demonstrate the important role played by the transverse gauge link in generating a 
gauge invariant answer. 

4.1 Feynman gauge 

First, we turn to the Feynman gauge. One-gluon exchange, see Fig. ^ trivially yields 

.(1) ^ (2p - k) ■ (2P - k) {2P - fc)+ 

' ((A;-p)2-m2 + z0)(fc2-A2 + 20) (A;+-iO)(P-A2 + zO) ■ ^ ' 

Employing the approximation (0), one finds 
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Figure 4: One-loop diagrams in the Feynman gauge. 



so that the small-x parton distribution reads to this order 



7r(x) = 2x J dp^p^ J ■ 



1 



(27r)2(27r)2(fc2 + A2-z0)2 \Z}(p-fc) D{p) j ' ^^^^ 
We now consider the one-loop contributions to the parton amplitude A2. In the small x region, 
the dominant contribution comes, as we already mentioned above, from diagrams with two-gluon 
intermediate states. There are six such Feynman diagrams shown in Fig. ^. As in Ref. 0, the 
gauge interaction between the D-meson and the rest is dominated by static Coulomb exchanges. 
It is easy to see that in the small-x limit, 

1 



A 



(2) 



ig 2p_- 



1 



where we have integrated over k2- and ^2 
interactions generate the amphtude 

. o 1 



(67) 



fk2_ 

(27r)2(fc? + A2-zO)(fc^ + A2-iO) 
On the other hand, the diagrams with two eikonal 



.(2) 
^2(c+d) 



ig 



d'^ko 



1 



k+-iOJ (27r)2 {kl + A2 - iO){kl + A2 - iO) ' 
Finally, we include the contribution for diagrams § (e) and (f) 

C2 1 [11 



.(2) 



-ig'^2p^ 



d^k. 



(27r)2(fc2 + ^2_^0)(fc^ + A2-z0) [/^(p-fci) D{p-k2 



(6J 



(69) 



where the minus sign indicates that the interference term is destructive. Summing over all con- 



tributions and using Eq. (62), we have the one-loop contribution to the parton amplitude, 



A 



(2) 



■9\ 
2! ^ 



rf2fc. 



1 



(27r)2 {kl + A2 - iO)(fc^ + A2 - zO) 

f 1 1 



D{p-k) D{p-k,] 
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D{p-k2) D{p) 



(70) 



(a) (b) 

Figure 5: One-gluon diagrams for the parton distribution in the hght-cone gauge. The dashed 
hne stands for absent hght-cone gauge hke and serves merely the purpose to guide the eye. 



where k = ki + k2 according to (|60|) . When the p ^ oo the amphtude vanishes as 1/p or as 
with r ^ when Fourier transformed. The argument is straightforward: The Wilson line 

represents an infinite-energy particle which has the opposite charge as compared to the spectator. 

When the distance between them goes to zero the gluon does not resolve this neutral composite 

system and the interaction vanishes. Eq. ([70|) agrees with the result of Ref. |0. 

Eqs. ( |65D and (^) are the first two terms in the expansion of the exponent. A computation of 

higher orders leads to the result for the amplitude 



^2(0) = iAp^J (fre'^ ^'^Ko [r^2xBP- + 

X j (fRS^-^'^-^^l-e^^2ias{KQ{{R + r)\)-KQ{R\)y^. (71) 



It has a factorized form of a convolution of the photon wave function Kq yry2xBP- + 'n^'^j and the 
scattering amplitude of the 00*-pair, separated by the distance r, of the struck quark, represented 
by the path-ordered exponential, and the spectator off the D-meson. In section ^, we reproduce 
this equation in the aligned-jet approximation for the dipole picture of '~f*D scattering. 

4.2 Light-cone gauge 

In this section, we calculate the same parton distribution in the light-cone gauge = 0. As we 
already emphasized in the introduction, this condition does not fix the gauge completely since the 
x_-independent gauge transformations leave it invariant. As a result the perturbation theory is 
not defined due to ambiguities in contour deformation of momentum integrals due to presence of 
an l/A;_,.-singularity in the gauge propagator 
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Most attempts in the literature to address the problem were pragmatic. A prescription was 
invented and then tested in practical loop computations. A justification by rigorous quantization 
procedures came later if ever ^. The fixing of residual gauge degrees of freedom, or maximal 
gauge fixing, can be done by means of boundary conditions on the gauge potential as is spelled 
out in the appendix 0. With this procedure one gets a set of regularization prescriptions on the 
light-cone pole corresponding to advanced, retarded and antisymmetric boundary conditions for 
the gluon potential, 

' Adv : 



[k, 



Ret: ^ . (73) 

PV- + 1 

Here only the pole associated with k^rif^ term, where the momentum k enters the vertex z/, of the 
gluon density matrix is displayed since only this part works in calculations performed later in this 
section. The other contribution can be obtained by hemiticity. 

In principle, any prescription for the light-cone singularity should work, so long as it is eas- 
ily generalizable to multi-loop calculations. Since all of the above prescriptions are not causal, 
loop integrations in particular Feynman diagrams can be pinched between poles and result into 
regularized spurious divergences. So that off-shell Green functions and renormalization constants 
do depend on the prescription at hand and, therefore, contain light-cone singularities. However, 



for a gauge invariant set of diagrams they cancel as they should, see, e.g., |Q for a two-loop 
computation with the principal value prescription. 

In our specific computation, the amplitude An, with the transverse gauge link, is gauge invari- 
ant by itself and independent on the regularization procedure. However, it is obvious that in the 
light-cone distribution (|T5|) the extra link does not contribute in the product of the parton am- 
plitude and its complex conjugate due to the unitarity condition [oo, oo; 0, cx5]"''[cxd, oc; 0, oo] = 1. 
In the following subsections, we verify these statements explicitly through the example of the 
previous section. 

To start with, let us address the tree level diagrams first. The contribution from the first term 
in the gluon propagator is the result of Fig. ^ (a) in the Feynman gauge. The term k^j^n^ yields 
a null contribution because fc^ contracted with the scalar D-meson vertex (P + P')^ vanishes by 
the on-shell condition. Contribution from Fig. ^ (a) is thus 

= TT, ^ . , I (2p " k) ■ (2P - k) - (2P - k)^k ■ {2p - k)- ^ ^ 



2(a) ((^_p)2_^2 + ^0)(A;2- A2 + i0) ^ ^ ^ ^ J+ [k+] \ ' 

(74) 

Next we exploit the simple fact 

k-{2p-k) = (p2 - m^) - ((A; - pf - m^) , (75) 
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Figure 6: One-loop diagrams in the light-cone gauge. 

where the first term dies out by the on-shell condition, in order to cancel the pion propagator in 
the amplitude with the second term. In the expression for the transverse gauge-link contribution 
we encounter a factor handled as 

p— ioofe+ 

= 2mx6{k+) , (76) 



with X depending on the prescription adopted for the pole 




X = < Ret : -1 . (77) 



As we expected, the advanced boundary condition, which sets the transverse link to unity, is 
supported by the perturbative calculation. Adding the contribution of the extra gauge link. Fig. 
1(b), 

= -#^A^2vr.x5(A:,) , (78) 

to (|7^, one gets the result of the Feynman gauge, Eq. (|6lD. 

Let us note in advance that in multi-gluon exchange diagrams, similarly to the one-gluon case, 
only kyU^ part of the gluon density matrix gives the leading contribution. The reason is that k^Uy 
part when contracted with the heavy meson line yields 




K-\2\P-Y:^k,yh]^={\P-Y.k,\ -M')-{\P-Y.k,\ -M'\, (79) 
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and since the dominant effect comes from the on-mass-shell D-meson intermediate states, both 
terms vanish independently. 



4.2.1 Advanced prescription 

Since the advanced prescription is implied by the advanced boundary condition A{oo) = 0, see 
the Appendix ^, the transverse gauge link at light-cone infinity does not contribute. This was 
observed in Ref. P] in a computation of a transverse momentum-dependent parton distribution. 



In this section, we perform an explicit calculation to reproduce the result of section in 
the light-cone gauge. Instead of showing the correspondence on a diagram-by-diagram basis in 
different gauges, we present the result of calculations for subsets of graphs. From diagrams (a) 
and (b) in Fig. ||, we obtain 

d^ko 1 



^(a+b) = ig^4:p_ J 



(27r)2(fc2 + A2-z0)(fc2 + A2-20) 

{D{p -k2)-D{p- k)) (Dip) - D{p - k^)) 



■ (80) 



D{p)D{p-k)D{p-k2) 
However, this is not the complete answer. In the present case, we have to take into account an 
additional diagram shown in Fig. |^. It is negligible in the Feynman gauge, but not in the light-cone 
gauge. The contribution of this graph is 

(fk2_ 1 D{p - fci) + D{p - fca) - D{p) - D{p - k) 

(27r)2(fc2 + A2-z0)(fc^ + A2-i0) D{p)D{p-k) 

(SI) 



The sum of the above two equations is exactly the right-hand side of Eq. (|7(]|) apart from terms 
which integrate to zero by the symmetry property of the integrand. 

The physics of the advanced prescription is that all final state interactions have been shifted 
to the initial state light-cone wave functions. Parton distributions are truly densities of partons. 
This scheme is the closest one to the original formulation of the Feynman parton model. 

4.2.2 Principal value prescription 

Now we discuss the principal value regularization for the light-cone singularity. One virtue of 
this prescription is that it does not introduce an absorptive part in the light-cone wave functions. 
Final state scattering phases then appear only through the extra gauge link derived in section 
|3.2| , as have been shown in the example of the single transverse-spin asymmetry in semi-inclusive 
deeply inelastic scattering 

A straightforward calculation gives the following expression for Fig. |^ (a) and (b), 

2^ f d'^k-i 



ig2p^l ^2ny{kl + \^-tO)ikl + x^-tO)\D{p-k) D(p-fc2)J ■ ^^^^ 
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Figure 7: One-loop diagram negligible in the Feynman gauge but contributing in the light-cone 
gauge. 



The diagram in Fig. |^ does not contribute. The transverse gauge link generates the diagrams (c), 
(d), (e) and (f) in Fig. H, whose straightforward computation gives, 

_ 2 f d'^k2 1 / 1 1 1 

^(c+d+e+f) - 2p. J J^2j^^^^—^j^j^^r^^—^^ - D{p-k,)l ■ ^^^^ 

Again the sum of the two contributions yields the same result as in Eq. ([TDD. 

The gauge link contribution in the parton distribution is unitary. When integrated over k, it 
does not contribute. Therefore, the square of the amplitude in Eq. (pSf ) must cancel against an 
interference between the two-loop gauge-link diagrams and the tree. This is indeed what has been 
observed in Ref. 0. The present discussion demystifies this cancellation. 

To conclude, the new gauge link found in section is necessary for reproducing the Feynman- 
gauge parton amplitude in the light-cone gauge calculation. However, it does not contribute to 
the integrated parton distribution. In particular, it does not play a role in modeling of leading 
twist "shadowing". 

4.2.3 Phase due to transverse link 

Actually, using the exponentiation property of soft radiation, one can compute the phase generated 
by the transverse link exactly in the limit of infinitely heavy mass of the D-meson, M —>■ oo. In this 
situation one can approximate the D-meson propagation by an eikonal line along its four-velocity 
= Pfj,/M. The parton amplitude reads 

^/lzl5° |d2j^e^(-P-^')-J^(p|T|[oot;,i?;-ooi;,i2][cx),oo;cx),0](/.(0,0)}|0). (84) 

Extracting the interaction of the 0-spectator with the target, we get 

J rf'i2e'(^--^')-^a(i?A)l^(i?A) . (85) 

where a stands for the (/)D-interaction as well as self-interactions, while 

ry/pxN _ {0\T{[oov,R;-oov,R][oo,oc;oo,0]}\0) 
vvy^A) (^Q\T{[oov,R;-oov,R]}\Q){0\T{[oo,oo;oo,0]}\0)' 
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represents the very effect of tlie transverse gauge link on the amplitude. In an abehan theory, the 
correlator of Wilson lines exponentiates exactly, 



W{RX) = e*'^(^^) , (87) 

where 



is given in terms of the gluon propagator (0|T {Afj_{x)A,,{y)} |0) = —iDn^{x—y) and the integration 
goes over the two pieces of path Ci = [oov, R; —oov, R] and C2 = [00, 00; cxd, 0]. A simple 
computation yields 

r d^k Q-ikR 

with 

In the propagator (^) only the term kyU^^ where k flows into the vertex of the transverse gauge 
link, works. Using Eq. ([77D, one finds 

(/^Adv(i2A) = , ¥^Pv(i?A) = -asMRX) . (91) 

This is an exact all-order result in the abelian theory, we are considering, and agrees with the lowest 
order predictions of Ref . . In a non-abelian theory, (f admits the web expansion [^, |27 



w 

where is the maximally non-abelian color factor and is a sum of corresponding web graphs. 
Eq. ( plD is the first term in the web expansion. 



5 Dipole picture 

This section is devoted to the demonstration of the equivalence of twist-two approximation for the 
structure functions at small values of the Bjorken variable and the aligned-jet kinematics of the 
dipole description of the 7*-hadron scattering. In the latter formalism, at small-xe, the virtual 
photon splits into a gg-pair long before the interaction with the target, due to large loffe time 
r ~ 1/(Mxb). On the second stage, the dipole interacts with the hadron. 

First we show this correspondence with the lowest order calculation in the scalar model alluded 
to above and then turn to an all-order proof in the semiclassical approach to small-xe QCD 
processes. 
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Figure 8: Dipole scattering off hadron target. 
5.1 Dipole scattering in scalar QED 

A simple calculation of the diagram in Fig. ^ within the scalar model discussed above produces 

d^k (fp Aal [1 1 1^ 



Ft^^'^x^, Q^) = J^' da aa{l - 2af J 



(92) 

where 

D{k) = aaQ^ + k'^ + - iO , (93) 

with a = 1 — a and momentum variables having the same meaning as in Eq. (|66D. Its Fourier 
transformation into the coordinate space gives 

P°''^(xB, g') = £ da ««(1 - 2a)2 J d^rK^ [r^aaQ^ + m^^ a{r) , (94) 

where the dipole cross section in two-gluon approximation reads 

a{r) = Aal J d^nl^Ko {{R + r)A) - Kq (RX) |^ + C(a^) . (95) 

Computation of the further terms leads to the exponentiation of the one-gluon exchanges similarly 
to Eqs. (p7D and (pSj). The dipole cross section can be expressed as a correlation function of two 
Wilson lines over the hadron states, 

J d^R (P'l (1 - U{R + r)U^{R)) \P) = 2n6{A+){2nf6^^\A)a{r) , (96) 

where A = P' — P. This makes the above observation on the exponentiation transparent. Here 
U is the abelian Wilson line 

U{x) = [00, x; —00, x] = exp ( —ig / dx-A^ (a;_, x) ) . (97) 



There are three regions in the dipole formula for the structure function (|94|) which lead to 
scaling. The small-r region, ~ ^/Q"^, corresponds to the genuine gluon contribution to the 
D-meson structure function. The aligned-jet regions ||28||: 



Q 
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is where the (anti)quark component dominates. So, to extract the leading twist-two quark contri- 
bution, we integrate in the vicinity of the end-point a = 0. Defining 



we get 




This result scales as a function of and thus defines the twist-two parton distribution, 

F^-2(xB, Q') = Ft''\xB, Q') . (101) 



5.2 Semiclassical approach 

In the previous section we have demonstrated in a model calculation that the dipole representation 
for the deeply inelastic structure functions in the aligned-jet kinematics is equivalent to the twist- 
two parton distribution functions at small momentum fractions. Here we generalize this result to 
all orders using the semiclassical picture of small-xe scattering in QCD. 

For the transversely polarized photon, the structure function has the form of a convolution of 



the photon wave function squared and dipole propagation though the target p9| , 

Ft^^'^xb, g') = 4iVc-^^ £ da aa (a^ + a^) J £r (r^aaQ^^ a{r) . (102) 

The dipole cross section, after extraction of total translations, cf. Eq. (PB|), is determined by the 
correlation of Wilson line, which represent the fast-moving quarks through the gluon field of the 
nucleon 



/ d'R ((^tr (l - U{R + r)U\R)))) , (103) 



with U being the non-abelian analogue of Eq. (|97|) with matrix valued gauge field and path ordering 
operation. At very high energies the color field is localized on the light cone with a profile in the 
transverse space. Namely, in the Feynman gauge, 

A^{z) = S{z.)Piz) , A_{z) = A{z) = , (104) 

it has the form of a shock wave [^, with the only non- vanishing components of the field- 



strength tensor G^±. A transformation to the light-cone gauge = with retarded boundary 
condition generates the pure-gauge color field of the McLerran-Venugopalan model [0, 

A{z) = e{z^)Vf3{z) , A4z) = 0. 
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Ft'''\xs,Q') = lQN,j^, j dp.p^ j SrKl[r^l2^Y{r). (105) 



In the aligned-jet kinematics (P5|), the structure function F2 = Ft + Fl can be written as 

i 

where we used the definition (^) and took into account that the longitudinal structure function 
is suppressed in this region, F^'^^^ ~ 

Now we are in a position to reproduce the same result ( |105[ ) computing the leading order and 
twist contribution to the same structure function, 

F*— 2(xb, g') = Xb (g(xB) - qi-XB)) , (106) 



at small xb, which is given in terms of the twist-two parton distribution 

Qix) = \l f^e»«-(P|^(0,0)[0,0;e-,0]7+^(e-,0)|P), (107) 

see Refs. [^, ^ |3^ for related considerations. Anticipating infinities due to unrestricted integra- 
tions, we introduce a regularized quark distribution 

q{x) = lJ^J dv-d'r]e'^^- {P'mv-,ri)[v-,r]-,^- + v-,rih+tp{^. + V-,ri)\P) , (108) 

which is related to the conventional one ( [107| ) by 

q{x) = 27r(5(A+) {2nf5^^^ ( A)g(x) , (109) 

where one sets P' = P in q{x). 

At very small x the bulk of contribution comes from the gluon component of the nucleon wave 
function while the constituent quark component is negligible. Therefore, we can approximate q{x) 
by computing the quark propagator in the external shock-wave field of the nucleon, 

q{x) = -lJ^J dil-d^ve'^^- (PlTr[r/_, ry; ^ + r^^, r7]7+z5(^ + v^,rj;r]^, ri)\P) , (110) 

where the trace is taken with respect to color and Dirac indices. The exact quark propagator in 
the shock- wave background ( |104| ) reads [0, |3^, |3^, |3l]]: 



iS{x; y) = iSo{x - y)d{x^yJ) (111) 
+ j d^z5{z,)tS,{x - ^)7- {9ix^)9i-y^)Uiz) - 9i-x^)9iy_)U\z)} tSoiz - y) , 
where 5*0 is a Fourier transform of the free quark propagator (pT]), 



So{x-y)= / ^e-^^-(^-^)So(fc) . 



Substitution of ( |111| ) into Eq. ( p,08| ) and subtraction of the free propagation yields 

g(x) = 16A^c^ j dp^p^ j d^rKl[r,j2^^ j d^R{P'\^ii (l - U{R + r)U\R))\P) . 

(112) 

After extraction of total translations on both sides with Eqs. (^6]) and ( |109| ), we reproduce the 
result of the aligned-jet kinematics ( [105| ). 
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6 DY versus DIS: universality of parton distributions 



In the Drell-Yan process, the gauge hnk in the parton distributions arises from the initial state 
interactions rather than from the final state. Using a procedure similar to the one carried out in 
section we can show that the transverse momentum-dependent parton distributions probed in 
DY production are, 

gDY(x, k) = y ^^e— «-+^'=-^(P|V^(^, ^)[-oo, oo; 7+[-oo, oo; 0, 0]c' ^(0, 0) |P) , 

(113) 

with the gauge link 

[-(X), oo;^_,^]c" = [-00,00; -cx),^][-cx),|;^_,^] . (114) 

The light-cone infinity ^_ = 00 has been replaced by ^_ = —00, reflecting the fact that the 
annihilation partner (antiquark or quark) from the other hadron starts its trajectory at ^_ = —00. 
As in deeply inelastic process we find again an extra contribution due to transverse components 



of the gauge potential which was not accounted for in recent treatments, see, e.g., Ref. |37 



In the light-cone gauge with the principal value prescription, the light-cone wave function 



is real. Therefore, the imaginary part in the parton distributions (51) and ( 113 ) arises from 
the interaction of the spectator with the gluonic field in the transverse link. Since in the DIS 
distribution the transverse link enters at ^_ = 00 while in DY at ^_ = —00 and the gauge 
potentials obey the antisymmetric boundary condition A(— 00) = — A(oo), we find that DY link 
differs by the sign in the phase from the DIS link, cf. Eq. (|48|), 



[— 00, 00; — 00, 0]pv = -Pexp ^— 2(7 J (i^-A(oo,^)^ . (115) 

Therefore, one can immediately understand in the lowest order calculation the change of the 
overall sign of the single transverse-spin asymmetry from DIS to DY case emphasized in Refs. 
TT| , |38[| . The phase of the transverse link, which is a genuine source of the imaginary part in 
the distribution for the principal value regularization of the light-cone gauge when is interacts 
with the target spectators, is opposite for DIS and DY transverse momentum-dependent parton 
distributions, see Eqs. (|48| ) and ( |115| ). This reflects a breakdown of the naive universality of 
unintegrated parton distributions: the Sivers function differs by an overall sign in DIS and DY 
p. 

When integrated over the transverse momentum, one flnds the following distribution from DIS, 
QDisix) = ^y'^e-^^«-(P|7/>(e-,0)[oo,0;e-,0]7+[oo,0;0,0]V'(0,0)|P), (116) 
and the distribution from Drell-Yan, 

QBYix) = ^J ^e-*^«- {Pm-, 0)[-oo, 0; ^, 0]7+[-oo, 0; 0, 0]^(0, 0)|P) , (117) 
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processes. The unitarity of the gauge Unk imphes that 



[(X),0;e-,0][oo,0;0,0] = [-00, 0; 0][-cx), 0; 0, 0] = [e-,0;0,0]. 



118) 



Therefore, both distributions are the same 39 
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|. In particular, this must be true 
in any model calculations. Therefore, although the approximation used in the derivation the DY 
scattering amplitude in Ref. |^ might be invalid for small k ~ xM, the integrated DY distribution 
in a full calculation must reproduce the distribution derived from DIS. 



7 Conclusion 

In the present study we have found that the transverse momentum-dependent parton distributions, 
as defined previously in the literature, are deficient in light-cone gauges which introduce additional 
singularities in perturbative calculations. We have derived their modification which consists of 
an extra transverse gauge link at the light-cone infinity. The regularization of spurious light-cone 
singularities can be related in certain cases to boundary conditions on the gauge field, some of 
which are consistent with the path integral quantization of the underlying gauge theory. These 
boundary conditions shift non-causal interactions, associated with nature of the gauge potential 
in the light-like gauge, into the initial or final state, or their mixture. One can get rid of the final 
state interactions in the light-cone gauge ^4+ = using the advanced boundary condition for the 
transverse components of the gauge field. With this choice, the extra gauge link in DIS parton 
distributions vanishes. Therefore, parton distributions are defined solely by the nucleon wave 
function and acquire the density interpretation of the naive parton model. Albeit, in this case the 
wave function is not real and apart from the structural information on the hadron it receives also 
the imaginary phase mimicking the final state interactions. The density interpretation is lost for 
other boundary conditions due to survival of the final state interactions. 

The application of these non-causal prescriptions might cause a problem in explicit compu- 
tations of multi-loop diagrams and one has to experience caution using them due to presence of 
non-integrable pinched singularities. The latter cancel however in gauge invariant quantities. 

Naive universality of transverse momentum-dependent parton distribution in deeply inelastic 
scattering and Drell-Yan lepton pair production is generally violated due to presence of residual 
gauge factors even in the light-cone gauge. Setting an advanced or retarded boundary condition, 
the transverse link in one of them vanishes, however, it persists in the other. However, the 
integrated parton distribution has to obey the universality condition due to unitary cancellation 
of gauge effects beyond the interquark separation points in cross sections. 
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A Maximally fixed light cone gauge 

Quantization of a gauge theory is hampered by an infinite redundancy of gauge degrees of freedom. 
To solve the problem one has to break the gauge invariance by imposing a condition on the gauge 
potential, = 0. Obvious requirements has to be fulfilled by the gauge fixing: (i) Every 

time one sets a constraint, one has to verify that there exists a choice of the gauge transformation 
matrix U, 

A^{x) ^ A'^ix) = Uix)A^{x)U\x) - -U{x)d^U\x) , (A.l) 

that it is realizable in practice, (ii) Next, one cannot gauge away a gauge covariant quantity like 
F^^ = d^Al - duAf^ - gf^'^A^^Al, so that -F[y4^] = must be consistent with this property. 

Since in the light-cone gauge one ascribes a specific (zero) value to the light-cone component 
of the gauge potential, it is very instructive to consider an equivalent example of temporal gauge 
fixing in a lattice gauge theory^ . Actually, one can consider merely a three-dimensional gauge 
theory on a cubic lattice, see Fig. ^. By a gauge transformation, one can eliminate one of the 
components of the gauge potential, say, 

A^{x,y,z) = 0. (A.2) 

This sets all link variables along the x-direction to unity, solid lines. This constraint does not fix 
the gauge completely since we are still free to perform x-independent transformations U{xo,y,z). 
To fix the residual freedom we can fix a subset of links oriented in the two-dimensional {y, z) 
plane. By the transformation t/(xo, y, z) we can fix 

Ay{xQ,y,z) = {], (A.3) 

but one cannot set at the same time Az{xo, y,z) = since otherwise we would set all link variables 
in the plane x = xq to unity and this would imply that we gauged away Fyz{xo) = 0. What we 
can do, however, is to fix with an (x, ?/)-independent matrix U{xQ,yQ, z) the component A^ on a 
line 

Az{xo,yo,z) = 0. (A.4) 

■^We would like to thank I.I. Balitsky for bringing it to our attention. 
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Figure 9: Maximal gauge fixing in three-dimensional lattice gauge theory. An example of the 
maximal tree Ax{x,y,z) = 0, Ay{xo,y,z) = 0, Az{xo,yo, z) = 0. Full links are set 



to 1 while dashed ones stand for conventional gauge links with non-zero gluon fields t/i,i+a 
exp {—igaAi). 



The supplementary conditions (|A.3|) and (|A.4| ) together with ( |A.2D provide the maximal gauge 



fixing for the gauge theory on a cubic lattice. Therefore, the maximal gauge fixing is determined 
by the maximal tree — a set of links such that adding one more results into a closed loop, 

forbidden by the gauge invar iance of the plaquette. 

In the continuum gauge theory the above residual conditions are translated into the boundary 
conditions on the gauge potential. Obviously, this is not the only choice for maximal gauge 
fixing. For instance, in the temporal gauge the auxiliary Coulomb fixing of residual gauge degrees 
of freedom d ■ A{to,x) = has been explored in Ref. |4^. It would be instructive to study a 



similar constraint for light-cone gauge A^{x) = where one may "freeze" A_[L_,x^,x) = and 
V ■ A{L_,x+,x) = 0. 

It is straightforward to extend the discussion of the lattice example to the light-cone gauge 
fixing in QCD. The maximal gauge fixing is achieved with a set of, e.g., retarded boundary 
conditions at L_ i = — oo: 

A+{x_,x^,x) = , A_{L_,x+,x) = , Ai{L_,0,x) = , A2(L_, 0, Li, X2) = . (A.5) 

The reason to fix the boundary condition at x+ = is that in the definition of the parton 
distributions all fields enter at x+ = 0. The maximally fixed light-cone gauge is obtained from 
other gauges by a series of gauge transformations. First, 

A+{x_,x+,x) = (A.6) 

with 

U\_{x) = [x_,x+,x; L_,x^,x] = P exp i^—ig J dx'_A^{x'_,x^,x)^ . (A. 7) 
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At light-cone boundary L_ one performs a gauge transformation on to nullify it, 

A_{L_,x+,x) = (A.8) 

with 

Ul{x) = [L^,x+,x; L^,0,x] = P exp l^—ig J dx'^A^{L^,x+,x)^ , (A. 9) 

etc. 

Exploiting Eqs. ( |A.5D , one can deduce the "inversion" formulas to express the gauge potential 
in terms of the field strength tensor, namely, 

dx _ [-^— ; '^-j- ; ^] ^ ; , X^ , , x"^ \^X , X ^ ; •^+7 ; ("A' -^^) 

+ y" 0, a;; L_, x^, a;]G_i(L_, a;)[L_, x'^, a;; L_, 0, as] , (A. 11) 

(ia;'„[L_, 0:4., a;; x+, a3]G'+2(a^'_, ic)[a;'_, x+, X] L_, a;+, a;] 
+ / (ia;',_[L_, 0, a;; L_, x^, a;]G'_2(-^-, a;)[L_, a;'_,_, a;; L_, 0, a;] 



+ / (ix'i[L_,0,Li,X2;L_,0,x'i,X2] 

xG'i2(L_, 0, x'l, X2)[L_, 0, x'l, X2; L_, 0, Li, X2] • (A.12) 

For x+ = 0, the transverse components simplify, 

Ai(x_,0,a;)= j dx'_G+i{x'_,0,x) , 

A2{x_,0,x)= dx'_G+2ix'_,0,x) + dx[Gi2iL_,0,x[,X2) , (A.13) 

and we dropped the path-ordered exponentials. Using these results we can compute the gauge 
propagator. For instance, 

{0\T{A4x.,x+,x)A{y_,0,y)}\0) = J^' dx'_ J'' dy'_{0\T {G+4x'_,x+,x)G+,{y'_,0,y)} \0) 

/d^k 
k2 + ,0 A;+ + 20 ' ^ ' 

where we have used the fact that the second term G12 in A2 does not work. Analogously, we get 
for other components 

d k ,1../^, „\ 1 ki 



{0\T{Mx.,0,x)A.{y^,y+,y)}\0) = t J — — e" 



(27r)4 fc2 + io A;+ - iO ' 

{0\T{A4x^,x^,x)A.{y^,y^,y)}\0) = ^ j 7^e-^'-^'~^)^^PV- , (A.15) 
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Figure 10: Feynman rules for the gluon propagator. 



where PV stands for the principal value prescription 

1 



PV±.i _ + 

k+ 2\k+ + i0 k 



(A. 16) 



This is actually the only components which contribute in practical calculation of section 4.2 



If we delinquently extend the results ( |A.14| ) and ( [A.15| ) out of the plane = 0, which 

implies that we exceed the amount of gauge fixing freedom in setting the boundary condition, we 
get the light-cone propagator 



\T{A,ix)A^iy)}\0) 



d'^k 



k^TLy 



kyTl^ 



(A. 17) 



(27r)4 A;2 + zOV"'"' k+ - k+ + iQ^ 

The gauge overfixing does not lead to complications in perturbation theory, since it boils down to 
the assumption G^y{x^ —oo,x^,x) 0. The direction of momentum flow k in Eq. (|A.17|) is 



assumed to be from the point x to y, see Fig. |T0|. The very same propagator was derived previously 
in the framework of path integral by means of a gauge transformation from the temporal gauge 
y4o = to the regularized gauge + ed+A^ = {e ±0) and canonical quantization in 
Ref. It was used more recently in the literature in the context of specific computations, see. 



e.g., Bli. 

If one uses the light cone gauge advanced boundary conditions L_ i = oo: 



A+{x-,x+,x) 



0, 



A^{L^,x+,x) = 0, Ai{L_,0,x) = 0, A2(L_, 0, Li, xs) = , (A.18) 



the contour deformation is moved into the opposite half-plane as compared to Eqs. ( A.14|) and 
( |A.15| ). With the same reasoning as preceded Eq. ( |A.17|) , one gets 



{0\T{A,{x)Ay{y)}\0) 



d^'k 



-ik-{y-x) _ 



1 



k^Tii, 



kufT'fi 



(A.19) 



(27r)4 k^ + ^oy^^^ k+ + iO k+ - iO ^ 

The conventional antisymmetric combination of boundary conditions ( [A.5|) and ( |A.18| ) bound- 
ary conditions, schematically, 

A(oo) + A(-oo) = (A.20) 

are not consistently implementable within the path integral approach. Moreover, one cannot 
find the corresponding gauge transformation which would result into this residual gauge fixing. 
However, if one accepts it in the "inversion" formulas, it results into the the propagator with l/k^ 
singularity regularized via the principal value prescription. 
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